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Abstract 


A stability margin evaluation method In terms 
of simultaneous gain and phase changes in all 
loops of a multi loop system is presented. A 
universal gain-phase ma. gin evaluation diagram is 
constructed by generalizing an existing method 
using matrix singular value properties. Using 
this diagram and computing the minimum singular 
value of the system return difference matrix over 
the operating frequency range, regions of 
guaranteed stability margins can be obtained. 
Singular values are computed for a wing flutter 
suppression and a drone lateral attitude control 
problem. The numerical results Indicate that this 
method predicts quite conservative stability 
margins. In the second example If the eigenvalue 
magnitude Is used Instead of the singular value, 
as a measure of nearness to singularity, more 
realistic stability margins are obtained. 

However, this relaxed measure generally cannot 
guarantee global stability. 
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Introduction 


Robustness criteria of multi-input 
multi-output (MIMO) feedback control systems using 
matrix singular value properties and their 
relation to the classical gain and phase margins 
have been a subject of research In recent 
years. I- 4 i n a multiloop system, the classical 
single loop ftyqulst tests are not adequate for 
gain or phase margin evaluation since the gain or 
phase is changed in one loop at a time. Using TFe 
minimum singular value of the system return 
difference matrix, criteria were developed in 
reference 4 for predicting guaranteed stability 
margins of multi loop systems in terms of either 
gain or phase change In all feedback loops. In 
this paper, the above method is generalized to 
include simultaneous gain and phase change in all 
loops. A universal gain-phase margin evaluation 
diagram is constructed. Using this diagram and 
computing the minimum singular value of the system 
return difference matrix over the operating 
frequency range, regions of guaranteed gain and 
phase margins can be obtained. Using this study 
and a singular value gradient expression derived 
herein, one can develop an algorithm for a direct 
design of robust multiloop control laws using the 
concept of cumulative constraint and multi- 
objective optimization. 5 - 6 

Singular values are computed for two examples 
for better understanding and evaluation of their 
relation to stability margins. The first example 
is a 29th order single-input single-output (SISO) 
system which represents an aeroelastic wind-tunnel 
wing model with a flutter suppression system. In 
reference 7 two flutter suppression control laws 
were synthesized that provided different degrees 
of robustness (gain and phase margins). These two 
control laws are examined using a classical 
Nyquist diagram and the singular value analysis. 
The correspondence between the stability margins 
predicted by the two methods are discussed. 
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The second example is an 8th order two-input 
two-output system which represents an experimental 
drone aircraft 8 with a lateral attitude control 
system. The phase and gain margins from classical 
single-loop tests and the singular value approach 
are determined and compared with the actual 
stability boundary. The plot of the return 
difference matrix eigenvalue in the complex plane 
is also studied. This plot provides useful 
qualitative gain and phase margin information not 
contained in singular values. 


Analysis 

Problem Formulation 

Let a multi loop feedback control system be 
described by the set of constant coefficient 
state-space equations (1) to (5). 


Tne stability can always be detected in the time 
domain by computing the eigenvalues of the closed 
loop system dynamic matrix 


F + G U LDH 

| G U LC ' 

BH 

r~ a 

i . 


and examining the sign of the real part of the 
eigenvalues. Traditional ly in single loop system 
design, the stability margins are evaluated as a 
gain or phase margin in the frequency domain using 
Nyquist diagrams. In a multiloop system, a 
frequency domain analysis may be formulated as 
follows. 

Write equations (1) to (5) and (8) in an 
augmented form as 


Plant 


x s = Fx s + G u u (1) 
z = Hx s (2) 
Control Law 

x c = Ax c + Bz (3) 
y = -[Cx c + Dz] (4) 
Interconnection 

u = r - y (5) 


Equation (1) represents an N s th order plant 
having Ng output measurements z modeled by 
equation (2) and N c control inputs u. 

Equations (3) and (4) represent an order 
feedback control law driven by the sensor 
output z. The reference input signal is r. In 
terms of a transfer function, the sensor output 
and the control laws are 

• z = CH (Is - F)-lG u ]u (6) 

u = [C(Is - A)-lB + D]z + r (7) 

respectively. Here s is the Laplace variable. 
Figure 1 shows the block diagram of a multiloop 
system with unity feedback in which a square 
diagonal matrix L is introduced at the plant 
input u to examine the system gain and phase 
margins. 

L = Diag[k n exp(j 4 , n )] (8) 


n = 1,2, . . .N c 

At nominal conditions, k n = 1 and 4> n = 0 
for all n and therefore L is the identity 
matrix. In a traditional single loop test, 
either kn or $ n is changed in one loop at 
a time. The gain or phase margin of the loop is 
defined as the smallest change in k n or $ n 
respectively from the nominal value for which the 
closed loop system remains stable. However, it 
would be more desirable to find the range in which 
both k n and $ n can be changed simulta- 
neously in all loops for which the system would 
remain stable. This is our present objective. 


x a = * a + BLu 
u = -Tx a + r 

where 
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H 
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( 10 ) 

( 11 ) 


( 12 ) 


Using equations (10) and (11), the plant 
input u and control output y can be written as 
functions of r in Laplace domain as 


u = (I '+ GfsjL)-^ 

(13) 

y = G(s)L( I + G(s)L)-lr 

(14) 

where G(s) = "C(Is - 7()-lTl 

(15) 

Using equations (6) and (7) one 
arrive at the relation (13) and (14) 

can also 
In which 

G(s) = -[C( Is - A)-lB + D] 


[H (Is - F)-lG u ] 

(16) 


Equations (15) and (16) can be shown to be 
identical. The transfer function matrix [I+G(s)] 
is defined as the system return difference matrix 
for nominal gain (i.e., L = I). For a steady unit 
sinusoidal input reference signal r, one can 
always compute [I+G(ju>)] over the range of 
operating frequency w. The multivariable system 
Is not robust to modeling errors if [I+G(ju)] is 
nearly singular at some critical frequency o> c . 

One can plot the determinant of [I+G(jw)3 In the 
complex plane and monitor its closeness to zero as 
it encircles the origin. But one cannot get a 
satisfactory notion of stability margin directly 
from the determinant locus because the closeness 
of a matrix to singularity cannot always be 
detected in terms of its determ1nant4 
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However, matrix singular values or principal 
gains are related to the system input-output 
magnitudes governed by equations (13) and (14) in 
the frequency domain, 2 » 3 and its properties are 
suitable to relate them to the stability 
margins. These properties are briefly reviewed 
next for later use. 

Review of Matrix Singular Value Properties 

Consider a sinusoidal input reference signal 
r = reJ ut for which the steady. state input u 
and output y are given by ueJ^ 1 _and _ 
yeJwt. The complex amplitudes u, y and r 
are related by the complex matrix relations 

u = T u (ju>)7 
y = T y (jio)r 

where T u = [I + G(ju)L]"' 

Ty = G(jw )L[I + GtjujL ]" 1 
= [I + (G(ju)L )- 1 ]- 1 


(17) 

(18) 

(19) 

(20) 


evaluation diagram as described in the following 
section. 

Universal Gain and Phase Margin Diagram 

Since the closed loop system is assumed to be 
stable for nominal gain (1*c«* L=I), the complex 
matrix T u must exist. This implies 

o(I+G) > 0 ( 29 ) 

The smallest singular value of the return 
difference matrix over the operating frequency 
range is a measure of stability margin. The 
objective of this section is to relate this to the 
maximum range of variation of k n and <j> n in 
the matrix L from their nominal values for which 
the perturbed system remains stable, i.e. 

o(l+GL) > 0 ( 3 °) 

It can be shown that the stability of the 
perturbed system is guaranteed if 

^(L- 1 -!) <o(I+G) (31) 


For convenience the argument (jw) will be omitted 
in the rest of the paper except when necessary. A 
scalar measure of amplitude of u is 

/ u*u = / r*[T u *T u Jr (21) 

where superscript * indicates complex conjugate 
transpose of a vector or matrix. The maximum and 
minimum singular values of T u provide a 
measure of the upper and lower bound of Ty in 
the sense 


assuming L“1 exists^. The basic reasoning 
is as follows. Since L“* and (I+G)’ 1 
exists, rewrite I+GL to separate L from 
the G matrix, as 

I+GL = [(L- 1 -I)(I+G)- 1 +I](I+G)L ( 32 

Since ( I+G) and L are assumed to be nonsingu- 
lar, (I+GL) is nonsingular if and only if 
[[(L-T-I) (I+G)“ 1 +I] is nonsingular. Condition 
(31) guarantees that 


max 


r*[T *T ]r 

'Vmax^V 


( 22 ) 


l[(L- 1 -I)(I+G)- 1 ]l 

max ^ 1 


or 


min 


r *W , , . „ (Tu) , 23 ) 


r*r 


The maximum and minimum singular values, denoted 
by 7 and o_, are given by the positive square 
root of the maximum and minimum eigenvalues 
of T U *T U (or T U V). The response . 
amplitude bounds hold true not only for sinusoidal 
Inputs but also for aperiodic and stochastic input 
signals. 2 Other important properties of 


singular values are 

o(G)a(L) > o(GL) ( 24 ) 

o(G) > |X(G)| > a(G) (25) 

If a(G) * 0 then G is singular (26) 

If exists 7(G) = 1/£(G"3) (27) 


If G is nonsingular, then a sufficient condition 
for (G + L) to remain nonsingular is 

7(L) <o(G) (28) 


Here G and L are any compatible square 
matrices and .X (G) denotes any eigenvalue of G. 
Using these properties and generalizing the 
analysis presented in reference 4, it Is possible 
to construct a universal gain-phase margin 


^[(L-i-lHI+G)- 1 ] < 1 ( 33 ) 

which ensures that [(L _1 -I)(I+G) _1 +I] is 
nonsingular. In arriving at equation (33), the 
inequality conditions (24), (28) and equation (27) 
are used. Note that this Is a conservative 
condition, and It is possible to construct 
an L matrix which violates this condition, yet 
fails to destabilize the system. 

If we consider simultaneous gain and phase 
changes in every loop using the L matrix given 
by equation (8), then 

c(L _1 -I) =yU-r) 2 + |-( 1 - cos 't'n ) {34) 

v *n n 


for all n with k n >0. Reference 4 considers 
only the classical cases of gain or phase changes 
keeping $ n B 0 or, k n =l, respectively. 

Using equation (34), the general case of 
simultaneous gain and phase changes can be 
examined. Equation (34) is presented In Fig. 2 as 
a universal diagram for gain and phase margin 
evaluation. For example if the smallest o_( I+G) 

for a system is 0.6 then the closed loop system 
will tolerate simultaneous gain and phase changes 
of — 1 . 5dB to +5.3dB, and -30° to +30°, 
respectively in all input loops. In a classical 
sense, when either gain or phase Is changed, the 
margins are -4.2dB and +8dB or +35°, respectively. 
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Using equation (34) the condition (31) can 
also be depicted as a plot of gain margin versus 
phase inaryin with smallest o_(I+G) as a 
parameter. This is shown in Fig. 3 to indicate 
regions of guaranteed stability. If the phase and 
gain changes k n and $ n are within the 
elliptic stability regions for a particular 
o(I+G), then the closed loop system must be 
stable. For example, if ojl+G) > 0.7 for all 
frequencies and 00° phase margins are required in 
all loops then guaranteed simultaneous gain 
margins are -2.6dB and +8.5dB, i.e., both the 
phase and gain can be changed in all loops in any 
manner within these limits, without destabilizing 
the closed loop system. 

From Figs. 2 and 3 it is also easy to verify 
the well-known result that for a Linear Quadratic 
Optimal state feedback problem the guaranteed gain 
margins are -6dB and + »dB and phase margins are 
+60° since o_(I+G)>l (see ref. 4). When both gain 
and phase changes are considered, for a given 
phase margin the corresponding gain margins can be 
established from Figs. 2 and 3. 

Computation of Singular Value and Derivative 

Since the frequency domain singular value 
analysis requires inversion of a large matrix 
(jlu-R) at a large number of frequencies (see 
equation (15)) an efficient computational 
method 9 * 18 is used. The basic idea is to 
transform "ft into an _upper Hessenberg matrix so 
that for all u> , (jloj-A) remains in upper 
Hessenberg form and the inversion problem 
(jIuj-E)“lB can be solved quickly by simple 
forward and backward substitution. Thus repeated 
upper and lower triangular transformations at 
each a) is avoided. 

The complex matrix (I+G) is of small order 
when multivariable systems involve only a few 
loops. It is possible to compute the singular 
value derivatives with respect to control law 
parameters with minimal additional computation and 
use them for direct design of robust multi loop 
control laws using the concept of multi objective 
optimization and cumulative constraints 8 * 8 . The 
expressions for singular value derivatives are 
similar to those of eigenvalue derivatives and are 
presented in appendix A. 

Examples 

Singular values are computed for two examples 
for better understanding and evaluation of their 
relation to stability margins. In the first 
example involving a single loop flutter 
suppression system, 7 the correspondence between 
the classical Nyquist diagram and the singular 
value plot are examined. Next a two input, two 
output system describing a drone aircraft 8 with 
a lateral attitude control system is analyzed in 
detail. 

29th Order Flutter Suppression System (SISO) 

In reference 7, 4th order flutter suppression 
control laws were synthesized for a 25th order 
state space system representing an aeroelastic 
wind-tunnel wing model. Two control laws were 
synthesized that provided different stability 
margins. The transfer functions for these two 
control laws are: 


Control Law (a) 

u = (-364.4) ( s-136.4) { s 2 +73-69s + 5697) deq 
z (s+2.057) (s+2057) ("s"^+’f6737s — + 2047) g 

Control Law (b) 

u = (1939. 4) (s+24.74) (s 2 +87. 63s + 13806) deg 

z rs+3.8^4) (s+3270)(S^+2b797"s“+ 1425) g 

Nyquist diagrams of the open loop transfer matrix 
G(s) are presented in Fig. 4. For control law (a) 
the gain margins are -4.1 dB and 2.6 dB and the 
phase margins are -22° and +41°. For control law 
(b) which is comparati vely robust, the gain 
margins are -5.0 dB and +12.3 dB, and the phase 
margins are -53° and +46°. 

A plot of the minimum singular value of (I+G) 
as a function of frequency with control law (a) 
and (b) are presented in Fig. 5. The points A, B, 
C and D on Fig. 5 correspond to the points on 
Fig. 4 where classical gain and phase margins are 
defined. If the singular values at point A, B, C 
and D are used with Mg. 2, the gain and phase 
margins are determined to be identical to those 
determined from Fig. 4. However, the minimum 
singular value occurs at a slightly different 
point denoted by E. If the singular value at 
point E (0.23 for control law (a) and 0.69 for 
control law (b)) is used with Fig. 2 the 
guaranteed gain or phase margins are determined to 
be -1.8 dB, +2.3 dB, or ±14° for control law (a) 
and -4.6 dB, +10.2 dB, or ±40° for control law (b) 
respecti vely. When point E is located in the 
Nyquist locus in Fig. 4, the worst direction of 
simultaneous gain and phase change are found to be 
+1.9 dB and -8.0° for control law (a) and -3.3 dB 
and +24° for control law (b) respectively. These 
values can be verified using Fig. 2. For this 
example the singular value plot of Fig. 5 
indicates only that control law (b) is more robust 
than control law (a). If the robustness is 
characterized by conventional gain and phase 
margins, quite conservative margins in certain 
gain or phase directions are obtained by the 
singular value approach. For realistic margins 
and information about the worst direction of gain 
and phase change, the Nyquist diagram is 
necessary. The Nyquist diagram in the SISO case 
is equivalent to plotting the eigenvalue of G 
versus frequency In the complex plane. Also 
X ( I+G) = 1+x ( G) . Thus qualitative gain and phase 
margin Information may be obtained for the MIMO 
case by plotting the complex eigenvalues of 
(I+G). This conjecture is examined in the next 
MIMO example. 

8th Order Drone Lateral Attitude Control System 

Tmtou 1 — 

Figure 6 shows the block diagram of the 
lateral attitude control system of a drone 
aircraft. 8 The plant state vector x s is 
defined as 


x s s L e ♦ $ ♦ $1 <$2j T 


The plant matrices, F, G, and H, as defined in 
Eqs. (1) and (2) are given in table 1. The 
control law matrices A, B, C, and D as defined by 
Eqs. (3) and (4) are given in table 2. The 
eigenvalues of plant x(F) are given in table 3. 
The eigenvalue at X = 0.1889 ± j 1.051 results 
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in an unstable Dutch roll mode. The elements of 
the input vector {uj ug} T are the elevon 
and rudder actuator servo commands respectively. 
All gain and phase changes are considered at the 
points denoted by X in figure 6. Figure 7 shows 
the Nyquist diagrams from classical loop breaking 
tests at the elevon loop and rudder loop 
respectively. From this diagram, the gain margin 
is 27.5 dB and the phase margin is -83.0° for the 
elevon loop designated as loop 1. For the rudder 
loop, the gain margins are -4.6 dB and •*« dB and 
the phase margins are +55° and -50.5°. This is 
the primary loop which stabilizes the Dutch roll 
mode and is designated as loop 2. 

Although the single loop values are often 
used as a measure of loop robustness, they may be 
inadequate to detect the overall system's weakness 
to simultaneous gain and phase changes in all 
loops or to unstructured uncertainties. To study 
this for the present two-loop problem, the minimum 
singular value of the return difference matrix 
(o(I+G)) over the operating frequency range Is 
plotted in Fig. 8. The minimum singular value is 
also the lower bound of the minimum eigenvalue 
magnitude |x_( I+G) I which is also plotted in 
Fig. 8. Th£ o_ is constant at 0.35 over low 
frequencies then drops to its lowest value of 0.25 
near 1.2 rad/sec which is close to the frequency 
of the unstable open loop pole. The o Increases 
sharply at higher frequencies and readies unity 
value asymptotically. Thus using the stability 
condition in Eq. (31), the stability_1s guaranteed 
if the L matrix has the property o(L _ I-I) 

<0.25. If the L matrix is structured as In Eq. 

8 which represents gain and phase changes in each 
loop, then by using Fig. 2, the guaranteed gain 
margins are found to be -2.0 dB and +2.5 dB. This 
means kj and k 2 can be changed from 0.8 to 
1.33 (with $1 =$2 a 0) in any manner without 
destabilizing the system. Similarly ^ and 
$2 can be changed between -15° and +15° (with 
ki=k 2 =l) without encountering instability. If 
both gain and phase are changed, one can again 
establish guaranteed gain margins for a given 
phase margin using Fig. 2 or 3. For example if 
the and $2 variation is between -10° 
and 10°, the kj and k 2 can be varied 
between 0.85 and 1.25 with guaranteed stability. 

The guaranteed stability margins and actual 
stability boundary for variations in kj and 
k 2 are compared in Fig. 9. The box ABDC 
Indicates the regions of guaranteed stability 
predicted from the minimum singular value. The 
curve IJ indicates the actual stability boundary 
obtained by computing the eigenvalues of matrix 
F a in Eq. 9 for many kj and k 2 and 
plotting the eigen-loci. The values of kj and 
k 2 where the real part of any eigenvalue Is 
zero determines one point on the IJ curve. Figure 
10 shows the same comparison when only $1 and 
4> 2 are varied. The actual stability boundary 
IJ is obtained by computing the determinant of 
[I+G(ju )] for many <|>i and $ 2 . On the 
boundary both the real and imaginary parts of the 
determinant are zero. The actual boundary can 
also be obtained by the eigen-loci method. Note 
in Fig. 10 that phase margin predictions from 
single loop f(yqu1st tests (Fig. 7) are Inadequate 
when simultaneous phase changes are considered. 
However, if the robustness Is characterized for 
design purposes by simultaneous gain or phase 
change tolerances at the input (Figs. 9 and 10) 
the singular value based gain or phase margin 


predictions are quite conservative, particularly 
in certain gain or phase change directions. 

In this example it is noted that considerable 
realistic but qualitative Information about gain 
and phase margins can be obtained from the eigen- 
values of (I+G) matrix, which is plotted in the 
complex plane in Fig. 11. Note that the minimum 
eigenvalue magnitude X(I+G)| plot in Fig. 8 
represents the radial dTstanci of the eigenvalue 
closest to the origin and is a measure of the 
closeness of (I+G) to singularity In a limited 
sense, i.e., for equal diagonal perturbation In 
the L matrix. If kj and k 2 are reduced 
simultaneously to 0.54, the eigen-diagram In Fig. 

11 shrinks radially about the point 1.0 without 
rotation and point C reaches the origin making the 
system unstable. Similarly both ^ and k 2 
can be increased indefinitely without encountering 
instability. If <}> j and $2 are changed 
simultaneously from zero to -38°, the diagram 
rotates about point 1.0 without distortion and 
point B reaches the origin making the system 
unstable. Likewise $1 and $0 can be 
increased simultaneously to +65* before point A 
reaches the origin. These specific stability 
boundary points can be verified from the actual 
stability boundary in Figs. 9 and 10. Although 
nothing can be said regarding stability from Fig. 

11 when gain or phase changes are unequal, it does 
provide qualitative information about the best and 
worst directions in gain or phase change. This 
useful design information cannot be obtained from 
singular value analysis. 

It is noted from Fig. 8 (as well as from 
Fig. 11) that the lowest Ixl Is 0.65 which 
occurs at a frequency of 1*54 rad/sec. Now if we 
are allowed to use the minimum eigenvalue magni- 
tude, A_(I+G) I, instead of the minimum singular 
value, o_(I+G), as a measure of nearness to singu- 
larity of the (I+G) matrix, then from Fig. 2 the 
relaxed gain margins are -4.4 dB and +9 dB. This 
means kj and k 2 can be changed In any 
manner between 0.6 and 2.8. Similarly from Fig. 

2, the relaxed phase margins are determined to be 
±38°. These two margins are plotted in Fig. 9 and 
10, respectively, and are denoted by the box 
EFHG. Comparison with the actual stability 
boundary IJ Indicates that these eigenvalue 
magnitude based predictions are more realistic. 
However, near the edge GH In Fig. 9 the relaxed 
boundary crosses over into the unstable region. 

The primary reason for this Is that the Inequality 
relation in Eq. (24) does not always hold true If 
singular values are replaced by eigenvalue 
magnitude. Hence, one cannot deduce Eq. (33) from 
Eq. (31) for eigenvalue magnitude except when all 
elements of a diagonal L matrix are equal. Thus 
the relaxed stability margin predictions using the 
condition 

jx(I+G)| 

generally cannot guarantee global stability. 


Conclusions 

A stability margin evaluation method for a 
multiloop system is presented by generalizing 
an existing procedure. The method Involves 
computing the singular values of the system return 
difference matrix over the operating frequency 
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range. The minimum singular value is related to 
the stability margins in terms of simultaneous 
gain and phase changes in all loops. Singular 
values are computed for two examples for better 
understanding and evaluation of their relation to 
stability margins. In the first example involving 
a single loop flutter suppression system of a 
wing, the correspondence between the classical 
Nyquist diagram and the singular value plot is 
examined. Next a two loop lateral attitude 
control system of a drone aircraft is examined in 
detail. Comparison with actual stability regions 
indicate that the minimum singular value based 
predictions are quite conservative. If the 
eigenvalue magnitude is used instead of the 
singular value as a measure of nearness to 
singularity more realistic stability margins are 
obtained for this example. However, this relaxed 
measure generally cannot guarantee global 
stability. The plot of the return difference 
matrix eigenvalues in the complex plane provides 
useful qualitative gain and phase margin 
information not contained in singular values. 


10. Peters, G. ; and Wilkinson, J. H. : Eigen- 

vectors of Real and Complex Matrices by LR and 
QR Triangularizations. Numerische Mathematik, 
Vol. 16, No. 3, 1970, pp. 181-204. 


Appendix A 

Consider a general complex matrix G whose 
singular values are a-j and the corresponding 
right and left normalized eigenvectors are vj 
and u-j , respectively. Hence by definition 


Gv-j = u-jo “j ( A1 ) 

G*uj = v-jo j 

i = 1, 2, ... N c (A2) 


The normalized eigenvectors satisfy the following 
orthogonal properties. 
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u i* u j = 6 i j vi*vj = j (A3) 

6 j j is the Kroneker delta which is unity when 
i=j and zero when i*j. Let p be a parameter for 
which derivative information is needed. Differen- 
tiating equation (Al) and (A2) with respect to p 
and then premultiplying the result by u^* and 
v-j*, respectively and adding them together, one 
obtains 


aG aG* av. 

u iVi + v i*ap~ u i + (U 1* G -V°1>3ir 


3U. 3a. 

+ (v.*G*-u i *a i )— = — (u^+v^v.) (A4) 


Using equations (Al) to (A3) in (A4), one obtains 


9q i = 1 (u,* 3£ v, + v.* ^G* u.) 

3P ? 3p 3p 1 

= real part of [u,* aG vJ (A5) 

3P 
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Table 1 Plant matrices F, G u and H for drone lateral attitude control system 


-.08527 

-46.86 

-0.4248 

0 

0 

0 


-0.0001423 

-2.757 

-0.06224 

1 

0 

0 


-0.9994 

0.3896 

-0.06714 

0 

0 

0 


0.04142 

0 

0 

0 

0 

0 


0.1862 

128.6 

-20.46 

0 

0 

- 20.0 


0 

0 

G = 0 


0 1 0 0 0 0 

0 0.07 1 0 0 0 


Table 2 Control law matrices A, B, C and D for 
drone lateral attitude control system 


0 0 
0 -2 


0.1491 0 

0 -4. 


D = 

0 

0 


0 

2.058 


Table 3 The eignevalues of the plant, x(F) 


-0.03701 

0.1889 ± j 1.051 
-3.25 
- 20.0 
- 20.0 


spiral mode 
Dutch roll (unstable) 
roll convergence 
elevon actuator 
rudder actuator 


L = 01AG 




| r -i f 

Tf+i- H(ls-F)' l G i^rT 

— 1 1 i L__ yj I I 

1 1 I i L 

i PLANT L 


-Cds-AI^B 
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PHASE MARGIN, + 0 .deg 
• — n 


MAXIMUM 1.2 
SINGULAR 
VALUE. 1,0 

OIL' 1 -!) 

0.8 



-4.0 -2.0 0 2.0 4.0 6.0 8.0 10.0 12.0 14.0 16.0 18.0 

GAIN MARGIN k, dB 



Fig. 1 Block diagram of a multi-input multi-output 
feedback control system with gain and phase 
change matrix at the input. 


Fig. 2 Universal diagram for multi loop gain-phase 
margin evaluation. 









GAIN MARGIN k dB 
n 

Fig. 3 Regions of guaranteed stability for various 
minimum singular value a_ (I+G). 



Fig. 6 Block diagram of a drone lateral attitude 
control system. 


phase mar gin 
♦ 41° AT 52 rad/sec 


PHASE MARGIN 
A + 44° AT 47 rad/sec 



Fig. 4 Nyquist diagrams of 29th order flutter 
suppression system from ref. 7 (arrows 
indicate increasing frequency). 


PHASE MARC'N 
A +55° AT 0.815 rad/sec 



Fig. 7 Nyquist diagrams from single loop breaking 
tests (drone lateral attitude control 
system). 



Fig. 5 Minimum singular value variation with fre- 
quency for 29th order flutter suppression 
system. 



Fig. 8 Variation of minimum singular value (a) and 
minimum eigenvalue (ixl) of (I+G) witF 
frequency (drone lat4ril attitude control 
system). 
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